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ABSTRACT. The use of molecular Rayleigh scattering for measurements of gas velocity and 
temperature is evaluated. Molecular scattering avoids problems associated with the seeding 
required by conventional laser anemometry and particle image velocimetry. The technique 
considered herein is based on the measurement of the spectrum of the scattered light. Planar 
imaging of Rayleigh scattering using a laser light sheet is evaluated for conditions at 30km 
altitude (typical hypersonic flow conditions). The Cramer-Rao lower bounds for velocity and 
temperature measurement uncertainties are calculated for an ideal optical spectrum analyzer and 
for a planar mirror Fabry-Perot interferometer used in a static, imaging mode. With this 
technique, a single image of the Rayleigh scattered light from clean flows can be analyzed to 
obtain temperature and one component of velocity. Experimental results sure presented for 
planar velocity measurements in a Mach 1.3 air jet. 


1. Introduction 

Laser diagnostics based on molecular scattering are receiving increased attention for 
use in aerospace test facilities. These methods can provide information directly 
related to the parameters of the molecular state of the gas including temperature, 
density, and bulk velocity. The elimination of the need for seeding is an important 
advantage of molecular scattering methods compared to conventional laser 
anemometry. Of the various molecular scattering techniques, Rayleigh scattering is 
the simplest and least expensive to implement. Rayleigh scattering has chiefly been 
used to obtain gas density by measuring the total Rayleigh scattered power. More 
information, however, is available from the spectrum of the Rayleigh scattered light. 

The Rayleigh scattering spectrum is directly related to the molecular 
velocity distribution function, which is a function of the thermodynamic parameters 
(temperature, density, pressure), composition, mean velocity, and turbulence 
parameters. Light scattered at a single scattering angle provides information 
related to one component of the molecular velocities. Thus one component of the 
mean velocity can be determined (from the shift of the spectral peak). The width of 
the spectrum is related to the gas temperature and turbulence. For low turbulence, 
gas temperature can be obtained; or, for a known temperature, turbulence intensity 
can be obtained. 

Extraction of the parameters of interest (temperature and velocity) from the 
spectrum can be accomplished using a sharp cutoff atomic or molecular absorption 
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filter [Shimizu et al.(1983); Miles et al.(1991)l, or with a high resolution 
spectrometer, such as a Fabry-Perot or Michelson interferometer. Point 
measurements of Rayleigh scattering spectra with scanning Fabry-Perot 
interferometers have been used to measure temperature in a hydrogen-air flame 
[Cattolica et al.(1976); Pitz et al.(1976)] and velocity and temperature in the 
exhaust of a small hydrogen-oxygen rocket [Seasholtz et al.(1991)]. Rayleigh 
scattering has been proposed for atmospheric temperature measurements [Schwiesow 
and Lading(1981); Rees et al.(1984)]. The application of single pass and multi pass 
scanning Fabry-Perot interferometers for high speed velocity point measurements 
has been previously analyzed [Seasholtz (1991)1. Imaging velocity measurements of 
seeded flows using a Michelson interferometer have also been reported by Seiler and 
Jager (1983). In this paper we examine an imaging technique to measure velocity 
and temperature using a planar mirror Fabry-Perot. The technique is an extension 
of ID temperatures and number density measurements in a heated air jet [Lock et 
al.(1992)] and is similar to Fabry-Perot imaging techniques for measurement of 
instantaneous velocity of solid surfaces [Gidon and Behar (1986)], temperature and 
velocity in the, upper atmosphere [Sivjee et al.(1980)], and solar corona 
measurements [Chandrasekhar et al.(1981)]. 

We first analyze a hypothetical planar imaging experiment where the 
measurement region is illuminated with a laser light sheet, and the Rayleigh 
scattered light is detected with a 2D array detector. It is assumed that the medium 
under study is a single component gas with no turbulence and contains no particles. 
The Cramer- Rao lower bounds for velocity and temperature measurement errors are 
calculated. We first calculate the error bounds for an ideal spectrum analyzer, 
which gives a basis with which to compare the performance of spectral 
measurements made with real instruments. The error bounds are then calculated 
for a planar mirror Fabry-Perot interferometer used in the static imaging mode. A 
specific example typical of hypersonic flow environments is discussed. 

The use of 2D Rayleigh imaging for velocity measurements is demonstrated 
with an experiment designed to measure the axial velocity component in a Mach 1.3 
supersonic free jet. 


2. Theory 

In this section we obtain theoretical lower bounds for uncertainty in velocity and 
temperature measurements for two cases where the scattered light is detected with a 
2D array detector. The first case assumes that an ideal instrument is used to 
measure the spectrum. The second case assumes that the spectrum is measured 
with a Fabry-Perot interferometer operated in a static, imaging mode (this is the 
technique described in the Experiment section). Only errors resulting from the 
Poisson statistics of the detected light are considered. Thus these results represent 
lower bounds for measurement errors. 

2.1. 2D RAYLEIGH IMAGING DIAGNOSTIC 

Consider a planar imaging experiment where a laser sheet uniformly illuminates a 
region in the object plane corresponding to the image of an array detector having N r 
rows (measured perpendicular to the beam propagation direction) and N c columns of 
pixels. Let the detector pixel size be L X *L X (square pixels) and let the system have a 
magnification M (i.e., the measurement plane has dimensions (N c L x xN r L x )/M). 
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The expected number of detected photons per pixel can be written [Eckbreth(1988)] 

< Npp ) = [hS] 

where Ei is the laser energy, A 0 is the laser wavelength, fdu/dfl) is the differential 
Rayleigh scattering cross section, n is the molecular number density, fi is the solid 
angle of the collected light, h is p lanck’s constant, c is the velocity of light, and e is 
the overall detection efficiency (including detector quantum efficiency and system 
losses). 

As an example, we consider air at 30km altitude (T = 230K, p=0.0118atm, 
n=3.83*10 23 m-3), A 0 =532nm, Ei=lJ, L x =23/rni, Nr=384, Mfi=0.05 (f/4 at M=l), 
e=0.05, and (d<7/dfl)=6.1*10' 32 m 2 /sr. For this case, equation 1 gives (N pp )=9.4 
detected photons per pixel. 

2.2. SPECTRUM OF RAYLEIGH SCATTERED LIGHT 


For a single component, low density gas with a Maxwellian velocity distribution, 
the normalized spectrum of the scattered light is [Fabelinskii (1968)] 

S(f) df = 2jg e-PMMJ-K-VP/a’K^ (2) 

where a = (2«T/m) 1/2 (m = molecular mass; k = Boltzman’s constant). The 
spectral peak is shifted from the laser frequency f 0 =c/A 0 by K*V/27T, where V is the 
mean velocity of the gas, and K=k s -ko (k s and k 0 being the wave vectors of the 
scattered and incident light). The magnitude of K is a function of the scattering 
angle 8 S and the wavelength A 

K= |K| = (47r/A)sin(0 s /2) (3) 

The observed Rayleigh scattering spectrum can be expressed as a set of counts 

<nj> =G,S(fj) Af (4) 

which are the expected number of counts in frequency interval fj to fj + Af; G R is 
the expected total number of counts. The spectrum (eq. 4) is thus a function of a 
set of three unknown parameters that can be expressed as the vector 

a = [G t ,V„Tj (5) 

where V K is the velocity component along K (i.e., V k =K-V/K). 

2.3. CRAMER-RAO LOWER BOUND 


The Cramer-Rao lower bound [Whalen(1971)J for the variance of the estimate of 
the parameter a-i is given by 

K«i) = Fin N (6) 

where T is the Fisher information matrix given (for p oisson statistics) by [Lading 
and Jensen(1980)] 


r 


ij - 


Y i a< n q ) n q > 

(n q ) da\ da j 


( 7 ) 


3 



If the parameters are uncorrelated, the estimate of the variance of a x is simply the 
inverse of the i th diagonal element of the Fisher matrix elements given by 


K«0 = 



( 8 ) 


2.3.1. Ideal Spectrum Analyzer. If we assume that the parameters are uncorrelated, 
we can easily evaluate the Fisher matrix elements and use equation 8 to get the 
Cramer-Rao bounds. Expressed as standard deviations a(ai)= [V(ai)] 1/2 , the lower 
bounds for measurement uncertainties for velocity and temperature are 

= » <K T ) = (2G^ji/2 ( 9 ) 

We now consider a subregion where it is assumed the parameters do not vary. The 
total number of photoelectron counts due to Rayleigh scattering detected from this 
subregion is then G E = (NPP) x number of pixels. The size of the subregion 
determines the spatial resolution of the measurements For example, assuming a 
20*20 array of pixels (which corresponds to. 0.5mmx0.5mm) and (NPP) =9.4 (as 
given in the above example), we have G E =3760 counts, T=230K, and a=370m/s. 
Thus the uncertainties in velocity and temperature are o(V K )=4m/s and <j(T)=5K. 
Note that the velocity uncertainty is proportional to the square root of temperature. 


2.3.2. Fabry-Perot Interferometer. The above error estimates were based on the use 
of an ideal instrument to measure the spectrum. In this section we perform the 
error estimation based on measurement of the intensity distribution of light imaged 
through a planar mirror Fabry-Perot interferometer. The expected number of 
detected photoelectrons (NDP) for the qth pixel is 


(NDPq) = (NPPq) J* j S(f)I t (f,0 r )df dApbel 


(10) 


where the Fabry-Perot transmission function [Vaughan(1989)] is 

W) = 1 + 'F 1 sin3(^/2) ( n ) 

with ip being the phase delay of the light for one pass through the interferometer 
given by 

^ _ Arydf cosflr 


and F = (2 Nr/7t) 2 , where N E is the reflective finesse. In this equation, f is the 
frequency of the light, /x is the refractive index of the medium in the Fabry-Perot 
cavity (taken as unity here), d is the Fabry-Perot mirror spacing, and 0 T is the 
angle between the ray and the optic axis. For simplicity, we neglected spectral 
broadening due to the finite size of the aperture of the collection optics. 

We now select, for this example, a particular Fabry-Perot configuration 
(mirror spacing d=15mm; finesse N E =20) and numerically evaluate the Fisher 
matrix elements 


rii - L 

q 



(13) 
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TABLE 1. Lower bounds for uncertainty in 
velocity and temperature measurements. 


Fringe 

number 

Fringe 

radius 

*(Vk), 

m/s 

*(T), 

K 

1 

38 pixels 

14 

44 

2 

69 

19 

49 

3 

90 

22 

61 

4 

106 

24 

71 

5 

121 

25 

68 

11 

185 

20 

81 

19 

286 

24 

83 

ideal 
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where the sum is over the same 20*20 pixel subregion used above. The results of 
this calculation are given in Table 1 for several locations of the subregion (denoted 
by fringe radii). The velocity uncertainty is about 20m/s (5 times that achievable 
with ideal spectrum analyzer) and the temperature uncertainty is about 40-80K (10 
times that of ideal). In practice, the assumption of constant finesse will not be 
valid; vignetting reduces the effective number of reflections in the interferometer 
cavity, thereby decreasing the finesse at larger fringe radii. This means that the 
finesse must be measured, as is done in the experiment described later in this paper. 

2.4. OTHER CONSIDERATIONS 

The example presented above considered an idealized measurement situation. There 
are always other factors that degrade measurement accuracy. For instance, if the 
gas is made up of an unknown mixture of component gases, the spectrum cannot be 
uniquely related to the temperature. If the flow is highly turbulent, and the 
measurement is obtained over a duration much greater than the turbulence time 
scale, the spectrum is broadened with both temperature and turbulence contributing 
to the spectral width. Thus, if one of these is known, the other can be determined 
from the spectral width, but both cannot be simultaneously determined. The 
Gaussian spectrum used in the example occurs only when the molecular 
mean-free-path is greater than the wavelength A=2x/K involved in the scattering 
process (i.e., y=A/mfp<<l). For values of y > 1, collective effects in the molecular 
velocities become important and must be taken into account in evaluating the 
spectrum shape. For y>>l, relatively simple continuum models can be used [Clark 
(1975)]; the spectrum in tlus case splits into three peaks and is referred to as 
Brillouin scattering. For the transition region, y~l, the spectral calculations are 
more involved and kinetic theory models, such as the Tenti S6 model [Tenti et 
al.(1974)] must be used. In general, ambient conditions fall into the transition 
region of y~l, while the lower densities typical of hypersonic flows have y<<l. 

Particles in the flow will strongly scatter. In some cases where the particle 
number density is low, the images can be processed to remove the obvious particles. 
If the particle number density is so large that the particles images overlap and 
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Figure 1. Experimental setup for measurement of velocity in Mach 1.3 free jet using imaging 
Fabry- Perot interferometer. 

cannot be removed, the spectral shape is changed from that of pure molecular 
scattering and temperature measurements cannot be obtained. It is important to 
note, however, that particles do not prevent the measurement of velocity (provided 
there is negligible particle lag) since the spectral peak due to particles is at the same 
frequency as the peak of the molecularly scattered light. 


3. Experiment to Map Velocity in Supersonic Free Jet 
3.1. EXPERIMENTAL SETUP 

To demonstrate the use of Rayleigh scattering for measuring velocity, we conducted 
the experiment shown in figure 1. The Mach 1.3 supersonic free jet (9.9mm exit 
diameter) was mounted on a linear translation stage so it could be moved along its 
axis. The beam from an argon-ion laser(2.5W, 488nm) was formed into a sheet 
(about 10mm high by 0.2mm thick) using a combination of cylindrical and spherical 
lenses. The light sheet crossed the flow at 45* . The receiving optics were aligned at 
45* to the flow direction and 90* to the propagation direction of the laser beam. 
Thus the K vector was aligned with the flow, and the axial velocity component was 
measured. A 250mm focal length f/3.6 lens was used to collect and collimate the 
scattered light. The collimated light was passed through a planar mirror 
Fabry-Perot interferometer (70mm diameter aperture and 14.96mm mirror spacing) 
and focused onto the CCD array with a 225mm focal length lens. The CCD array 
had 384 rows by 576 columns of 23pm x 23pm square pixels (giving a field size in 
the flow of 9.8mm* 14.7mm) and was cooled to about -50* C to reduce dark current. 
The CCD pixel data was digitized with 12 bit resolution and transferred to an 80386 
PC via an IEEE-488 interface for storage and analysis. 

The initial step in the measurement process was to measure the instrument 
function [It(^), for f=f 0 J. A diffusely scattering target was placed in the field of 
view and illuminated with laser fight scattered from a rotating diffuse target (the 
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Figure 2. Pattern of subregions used for data analysis (5x5 pixel subregions shown; actual size of 
regions used was 20*20 pixels). 

rotating target was used to reduce speckle in the recorded image. The Rayleigh 
scattered light from the flow was then recorded. Exposure times were 0.1s for the 
instrument function and 1.0s for the Rayleigh scattering. 

Although the air supplied to the nozzle was filtered, it was not feasible to 
supply filtered air to the air surrounding the jet. Particles in the ambient air were 
thus mixed into the flow, and temperature measurements were not realizable (even 
though the technique described by Kay and Sadler (1991) was used to reduce the 
effect of particles). 

3.2. DATA REDUCTION 

The data reduction was done in two steps. First, the instrument function image was 
analyzed to determine the center of the circular fringe pattern, the finesse, and the 
phase corresponding to the laser frequency. Because the effective finesse and phase 
varied over the image plane, it was necessary to determine them as a function of 
location in the image (typically, the finesse varied from about 20 at the center to 
about 4 at the edges of the image). About 500 20*20 pixel subregions were used. 
An integral number were positioned around the each of the inner fringes and around 
every other outer fringe (so the subregions did not overlap). A typical pattern of 
subregions is shown in figure 2. A 4“* order, 2D polynomial (15 parameters) was 
then fit to these values of phase and finesse. The coefficients for the polynomial fits 
were saved for the second step (analyzing the image of the Rayleigh scattered light). 

In the second step, the image of the Rayleigh scattered light in each 
subregion was fit to a model function given by equation 10 using the S6 model of 
Tenti et al.(1974) with the gas properties of air. The velocity data obtained at each 
subregion were fit to a 2D 4 th order polynomial, which was used to generate contour 
plots of constant axial velocity. 



3.3. RESULTS 


Figure 3 shows constant axial velocity contours for 5 axial locations along the jet 
axis (X=20,30,50,70, and 90 mm). Recall that these planes are oriented at 45* to 
the jet axis. At the locations closer to the nozzle exit, the maximum velocity is in 
good agreement with the 423m/s velocity calculated from the measured total 
pressure and temperature (assuming isentropic flow). As the measurement plane 
moves downstream, the maximum velocity decreases as the core flow mixes with 
ambient air. The contours are elliptically shaped because of the orientation of the 
measurement plane. We also can see that the center of the jet is shifting to the 
right as it moves downstream. This is probably caused by a slight misalignment of 
the jet axis. 

The uncertainty in the velocity estimates due to photon arrival rate statistics 
is quite small, less than lOm/s. Errors due other sources most certainly outweigh 
these statistical errors. One potential error source is the drift in the Fabry-Perot 
mirror spacing between recording the instrument function and the Rayleigh data. 
The r.m.s. deviation of the contours Rom the individual subregions data was about 
10-15m/s. 

The data reduction requires rather lengthy processing, chiefly because of the 
numerical integration in equation 10. The time per iteration for each subregion is 
about 15s using a 80386/33MHz/Weitek 3167. Thus an image with 500 subregions 
and an average of 5 iterations, requires about lOhr to process. No attempt was 
made to optimize the code at this time, so it should be possible to substantially 
reduce the data reduction time. 


4. Concluding Remarks 

The results of this work show the feasibility of using measurements of the spectrum 
of molecular Rayleigh scattering as a diagnostic suitable for many high speed flow 
experiments. The technique is capable, under certain restrictions, of mapping gas 
density, velocity, and temperature. 

The experimental work reported here used a CW argon-ion laser, which 
results in time-averaged measurements. Use of a pulsed laser and gated detection 
would offer several advantages. Measurements based on a image from a single laser 
pulse would give 2D maps of the gas parameters at an instant of time. The effect of 
background light, such as flame luminosity, could be greatly reduced. And, particle 
images, if not too dense, should be easier to identify and remove from the flow field 
image. 

A large number density of particles in the flow generally will prevent 
temperature and density measurements. It is important to note, however, that 
velocity can be obtained with or without particles, provided the particles follow the 
flow. Care should be taken to prevent spurious laser light from reaching the 
receiving optics. Even though this light can be accounted for in the data reduction, 
it will degrade the signal-to-noise ratio and measurement accuracy. This problem 
is particularly severe in internal flows where it is very difficult to eliminate all 
spuriously scattered light. The technique described by Miles et al. (1991) for 
blocking light at the laser frequency should be useful for high velocity flows where 
the frequency of the Rayleigh scattered light is shifted well away from the laser 
frequency. 
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Figure 3. Measured constant axial velocity contours. Distances from nozzle exit are 
X=20, 30,50, 70, 90mm. Isentropic velocity = 422m/s. 
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